Non-linear unbalanced Bessel beams: 
Stationary conical waves supported by nonlinear losses 
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Nonlinear losses accompanying self-focusing substantially impacts the dynamic balance of diffrac- 
tion and nonlinearity, permitting the existence of localized and stationary solutions of the 2D+1 
nonlinear Schrodinger equation which are stable against radial collapse. These are featured by lin- 
ear, conical tails that continually refill the nonlinear, central spot. An experiment shows that the 
discovered solution behaves as strong attractor for the self-focusing dynamics in Kerr media. 

PACS numbers: 42.65.Re, 42.65.Tg 
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One of the main goals of modern nonlinear wave 
physics is the achievement of wave localization, station- 
arity and stability. While in a one-dimensional geometry 
(e.g., in optical fibers), nonlinearity suitably balances lin- 
ear wave dispersion, leading to the soliton regime, in the 
multidimensional case, nonlinearity drives waves either to 
collapse or instability. In self-focusing of optical beams, 
for instance, many stabilizing mechanisms, such as Kerr 
saturation, plasma-induced defocusing, or stimulated Ra- 
man scattering, have been explored, and are being the 
subject of intense debate, mainly in the context of light 
filamentation in air or condensed matter Q . These mech- 
anisms, however, are either intrinsically lossy, or, due to 
the huge intensities involved, are accompanied by losses, 
which lead ultimately to the termination of any soliton 
regime. Similar pictures can be traced in all phenom- 
ena commonly discussed in the context of the nonlinear 
Schrodinger equation (NLSE), as Bose-Einstein conden- 
sates (BEC) or Langmuir waves in plasma 0- Nonlinear 
losses (NLL) arise in BEC from two- and three-body in- 
elastic recombination, and as the natural mechanism for 
energy dissipation in Langmuir turbulence 0. 

The question then arises of whether any stationary and 
localized (SL) wave propagation is possible in the pres- 
ence of NLL. The response, as shown in this Letter, is 
affirmative. These SL waves cannot be ascribed to the 
class of solitary waves, but are instead nonlinear conical 
waves (as the non-linear X waves 0]) of dissipative type, 
whose stationarity is sustained by a continuous refilling 
of the nonlinearly absorbed central spot with the energy 
supplied by linear, conical tails. These waves are not only 
robust against NLL, but find their stabilizing mechanism 
against perturbations in NLL themselves. 

Among the linear conical waves [f| , the simplest one is 
the monochromatic Bessel beam (BB) 6], made of a su- 
perposition of plane waves whose wave vectors are evenly 
distributed over the surface of a cone, resulting in a non- 
diffracting transversal Bessel profile. Despite the ideal 
nature of BBs (they carry infinite power) , they not only 
have revealed to be a paradigm for understanding wave 



phenomena, but also have found applications as diverse 
as in frequency conversion, or in atom trapping and align- 
ment 7]. Of particular interest for us is the finding Q 
that the BB is describable in terms of the interference of 
two conical Hankel beams |5j , carrying equal amounts of 
energy towards and outwards the beam axis, and yielding 
no net transversal energy flux in the BB. 

What we demonstrate here is that a superposition of 
inward and outward Hankel beams with unequal am- 
plitudes, i.e., an "unbalanced" Bessel beam (UBB), de- 
scribes the only possible asymptotic form of the SL, non- 
singular solutions of the 2D+1 NLSE with NLL. For 
this reason, we call the solution "non-linear UBB" (NL- 
UBB). We then show that NLL-UBB solutions do exist 
and are stable against radial perturbations, in the impor- 
tant case of Kerr nonlinearity with NLL. Unbalancing, 
that manifests as a reduced visibility of the radial Bessel 
oscillations, creates the required inward radial energy flux 
from the conical tails of the beam to refill the nonlinearly 
absorbed central spot, whose transversal pattern depends 
on the specific nonlinear phase (nondissipative) terms in- 
cluded in the NLSE. Contrary to linear conical waves, 
achievement of stationarity by refilling imposes generally 
a lower bound to the cone angle of the UBB. 

In a self-focusing experiment in water, we demonstrate 
the self-generation of a NL-UBB from a Gaussian wave- 
packet, which evidences that the NL-UBB acts as a 
strong attractor in the self-focusing dynamics, and hence 
a significant role of NLL in light filamentation. 

To start with, we consider the 2D+1 NLSE 
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for the propagation of a light beam E — Aexp(— iiot + 
ikz) of frequency w in a Kerr medium (other nonlinear 
phase terms could be included as well) with NLL. In 
Vj_ = (d x ,d y ), k = nuj/c is the propagation con- 
stant (with n the refraction index and c the speed of 
light in vacuum), n-i is the nonlinear refraction index, 
and 0^> > (K = 2, 3, . . .) is the multi-photon ab- 
sorption coefficient. For the real amplitude and phase 



2 



[A = aexp(itp), a > 0], yields 
1. 
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Stationarity of the intensity profile (d z a 2 = 0) requires, 
from J2J|, if = 4>{x, y) + g{z). Then J2J imposes the linear 
dependence g(z) — —dz, where 5 is a constant wave vec- 
tor shift. Equations J2J and © then lead, for cylindrical 
beams, to the eigenvalue problem 
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[prime signs stand for d/dr, with r = (x 2 + y 2 ) 1 ^ 2 ], with 
boundary conditions a(0) = a > 0, a'(0) = 0, 4>'(0) = 
0, and the requirement that a(r) — > as r — > oo for 
localization. Eq. © establishes that in a SL beam, the 
power (per unit propagation length) entering into a disk 
of radius r must equal the power lost N r within it. In 
absence of NLL (@( K ' = 0), this condition demands plane 
phase fronts ((/>' = 0) and no radial energy flux. This is 
the case of the sech-type Townes profile in Kerr media 
0, associated to wave vector shift 6 < 0, and of we akly 
localized Bessel beams in linear |(j or Kerr media |l0j . 
with infinite power and wave vector shift <5 > 0. 

With NLL, instead, the conditions of refilling JSJ and 
localization [a(r) — > as r — > oo] require an inward radial 
power [lhs of |J5J|] that monotonically increases [rhs of JSJl] 
with r up to reach, at infinity, a constant value equal to 
the total NLL, iV^, assumed they are finite. Stationarity 
with NLL is thus supported by the continuous refilling of 
the more strongly absorbed inner part of the beam with 
the energy coming from its outer part. Phase fronts can- 



not be plane, since 



-kNoc/2-Kra as r — + oo. As for 



the amplitude, the asymptotic value of 4>' and the change 
a(r) = b(r) / y/r in |@}, lead, when retaining only the slow- 
est decaying contributions, to the Newton-like equation 
b" = ^kSb + PN^/An 2 ^, that represents the "motion" 
of a particle in the potential V(b) = kSb 2 + k 2 / %ir 2 b 2 . 
Since bounded trajectories b(r) [leading then to a — > 0] 
under this potential can exist only for strictly positive 5, 
we conclude that SL waves in media with NLL can only 
have positive wave vector shift, S > 0. The solution of 
the Newton equation then yields the asymptotic behav- 
ior a (r) = {{c 1 +c 2 cos(2\/2k5r + c 3 )]/r} 1 / 2 , with c x > 0, 
c i — c 2 = ^^L/8tt 2 (5, that represents radial oscillations 
of contrast C — \c 2 \jc\ about an equilibrium point that 
approaches zero as 1/ y/r. SL beams in media with NLL 
carry then infinite power, and have superluminal phase 
velocity (5 > 0). These asymptotic features are more 
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FIG. 1: (a) For pure NLL (K = 4), radial profiles a 2 /, 
SL beams with decreasing 8/f3 {K ) I*' 1 = oo, 0.33, 0.25, 0.20 
(from lower to higher ones), and their asymptotic UBB. Nor- 
malized radial coordinate is p — \p2Mir (b) Allowed cone 
angles 8 = y/26/k versus peak intensity Jo normalized to 
the intensity (fc//3 (x) ) 1/(K ~ 1) in the cases of pure NLL, and 
NLL+Kerr (numerically calculated and from Eq. [0J|. (c) For 
NLL+Kerr, amplitude profile a/a for S/P^I^ 1 = 0.24 
and knzIo/nS = 6.07, and its asymptotic UBB. The pure Kerr 
case is also shown, (d) Eigenvalue spectrum of the perturba- 
tions with m = to the SL solutions with Kerr nonlinearity 
g = 1 and NLL (K = 4) 7 = 0, 0.1, 0.2. 



meaningfully expressed in terms of the UBB 
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formed by two nondiffracting Hankel beams of the first 
and second kind 5], of same cone angle 9 = y/26/k, but 
different weights, a ou t and aj n , that must be related by 



ao(Kn| 2 - KutHA = N * 



(7) 



The balanced superposition of Hankel beams (e.g., 
a out = a in = 1) J us t gives the original, nondiffracting 
Bessel beam aoJo(V2kS r) exp(-iSz), with no net radial 
energy flux, and with oscillations of maximum contrast 
(7=1. Instead, unbalancing creates an inward radial 
power that is manifested in a lowering of the contrast 
(7= M/ci = 2|a in ||a ou t|/(|a in | 2 + |a out | 2 ) of the Bessel 
oscillations, reaching (7 = (no oscillations) in a pure 
Hankel beam. 

We stress that this analysis holds irrespective of the 
nonlinear phase terms in the NLSE (Kerr nonlinearity, 
Kerr saturation, . . . ), since they rely on the only assump- 
tion of finite NLL. We can therefore state that the conical 
UBB represents the only possible asymptotic form of SL 
waves in nonlinear media when the effects of NLL are 
taken into consideration. Note that the case of linear 
losses is excluded, since Noo = 00 for K = 1. It should 
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be also clear that the actual existence of a SL solution 
of the NLSE, and the characteristics of its linear asymp- 
totic UBB (cone angle, a- m and a ut); depends on the 
particular nonlinear phase terms in the NLSE. 

We first solved numerically (@J and JHJ without the 
Kerr term, to appreciate the effects of pure NLL [Fig. 
ma)], an d found that SL solutions [a(r) — ► 0] exist indeed 
with any peak intensity Iq 



Oq and wave vector shifts 
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where g K = 1.67, 0.27, 0.19, 0.16 . . . for K = 2, 3 . . . 
As seen in Fig. m & )j shortly away from the central peak, 
the radial profile becomes undistinguishable from that 
of the UBB of same 5, matched NLL Noo and contrast 
C . Moreover, pure NLL creates a gap in the allowed UBB 
cone angles 9 = y/26/k that increases with intensity [Fig. 
mb)], and that becomes significant at intensities com- 
parable to the characteristic intensity (k/ft^') 1 ' ^ K ~ 1 >. 
Intuitively, this gap originates from the fact that the SL 
profiles widen and delocalize [Fig. m a )] as the cone angle 
diminishes. This causes the total NLL, Noo, to increase, 
a situation that cannot be sustained down to the limit of 
zero cone angle (5 — > 0), which would not allow for any 
radial energy flux. 

The extension to include nonlinear phase terms in the 
NLSE can be readily understood from the case of pure 
NLL. Figure m^O shows the permitted cones angles in 
the case of focusing Kerr nonlinearity (712 > 0) , obtained 
from numerical integration of (@J and j£| with K = 4 
(similar results hold for other values of K). The modi- 
fication of the allowed region of cone angles can be at- 
tributed to the nonlinear phase shift at the central spot. 
In fact, assuming that the existence of a localized solu- 
tion is now determined by the effective wave vector shift 
^cff = ^ + ^nb with 8 n \ = kn 2 Lo/n for Kerr nonlinearity, 
we replace 6 with <5 e g- in Eq. |JSj, to obtain 



6 > maxj^/^/^- 1 - kn 2 I /n, oj 
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(where we set to zero negative values) as an accurate 
expression for the allowed linear wave vector shifts in 
Kerr media [see Fig. m°) for the associated cone angles 
9 = ytW/k]. Similar expressions as © can be obtained 
for other nonlinear phase effects. 

When NLL dominate over Kerr nonlinearity [right part 
of Fig. or 9Kf3 {K) Io^ 1 > kn 2 Io/n 2 } the SL pro- 

files (not shown) do not substantially differ from the case 
of pure NLL. For the Kerr-dominated case [left part of 
Fig. mb), or gK^L*- 1 « kn 2 L /n 2 }, Fig. [flc) shows 
a representative SL profile. The central peak and in- 
ner rings are nearly identical to the Kerr-compressed, 
Bessel-like beam in lossless Kerr media 0]j as seen in 
Fig. m c ); though the small NLL leads to a slight con- 
trast reduction. The central peak can be approached 
by the Bessel beam a Jo(^/2fc<5 e fjr) and hence its 



width by 1/ ^2k8 e g. Outer rings gradually shrink up to 
become phased with the asymptotic UBB of wave vector 
shift 5 [Fig. m c )]- Note that, as opposed to the pure 
NLL regime, the beam does not widen indefinitely as the 
cone angle diminishes down to its lower bound (9 — ► + 
in the Kerr-dominated case), but is limited to a maxi- 
mum beam width 1 / y j / 2kS e ff = y // n/2k 2 n 2 Io. This fact 
entails, contrary to the pure NLL regime, a limitation to 
Noo as 9 — > + , explaining why the unbalance mechanism 
for replenishment can support stationarity at arbitrarily 
small, but positive cone angles. 

We also studied the stability of the NL-UBB solutions 
against perturbations. With the dimensionless quantities 
p = \/2k5r, £ = Sz and A — A/ao, we rewrite Q as 
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where X7 2 = d 2 + (l/p)d p , g — uin 2 al/cS and 
7 = (3 (K ' > al K ~ 2 /2S. Following a standard Bogoliubov- 
deGennes procedure 0] > we introduce a perturbed solu- 
tion 
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where Ao(p, £) = a{p) exp[z</<(p)] exp(— z£) is a SL solu- 
tion of (|10fl . 9 is the polar angle, and m = 0, 1 . . ., into 
(|10f) . to obtain, upon linearization, the (non-self-adjoint) 
eigenvalue problem 
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where H = — V 2 + m 2 / p 2 — 1 — 2ga 2 , with boundary con- 
ditions u(p) —>■ 0, v(p) — * as p — > 00. The existence of 
an eigenvalue SI with Imfi > would render unstable the 
SL solution A®. To obtain numerically the eigenvalues, 
we transform l|12H into an algebraical eigenvalue prob- 
lem by introducing a mesh of small size s in the range 
[p = 0, p = L] with large L, writing the differential op- 
erators as finite differences, and imposing the boundary 
conditions u(L) = v(L) = 0. The numerical diagonal- 
ization of the 2N x 2N matrix (jV being the number of 
mesh points) of the eigenvalue problem provides a set of 
2N eigenvalues. The behavior of the continuous distribu- 
tion of eigenvalues of (|12|l is inferred by extrapolating the 
results to the limit s — > and L — > 00. We investigated 
increasing N up to 4000 (L = 400, s = 0.1), as limited 
by the accessible memory of our computational facility. 

Figure m^) shows a typical eigenvalue spectrum for 
radial perturbations (m = 0) to three SL solutions with 
same Kerr nonlinearity (g — 1) and increasing NLL. 
Eigenvalues with ImSl = are not shown for clar- 
ity. Clearly, the effect oj NLL is to decrease the imag- 
inary part oj the eigenvalues, driving the system towards 
stability: The Bessel-like solutions in pure Kerr media 
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FIG. 2: Measured fluence profile (dots), fitted UBB, and half 
Bessel (left) and Airy (right) profiles. 



(squares) are unstable. At NLL strength 7 = 0.1 (stars), 
the positive imaginary parts are strongly reduced, and at 
7 = 0.2 (circles) no signs of instability are present. We 
performed the same analysis for m = 1, 2, . . . For dipolar 
perturbations (m = 1) no instability emerges neither in 
the pure Kerr nor in the Kerr + NLL model. In contrast, 
both systems turned out to be unstable for quadrupolar 
and higher-order perturbations, a result that can be re- 
lated to the fact these perturbations involve modulation 
far from the central spot, where intensity is weak and so 
NLL cannot play its stabilizing role. 

The experiment that we present demonstrates that the 
NL-UBB stationary solution indeed acts as strong attrac- 
tor for the transient dynamics of light beam self-focusing 
in (weakly dispersive) Kerr media. Recently |12( we have 
shown that light filaments in water do not behave as 
soliton-like beams, but they spread after being clipped 
by an aperture and reconstruct themselves after being 
blocked by a stopper, as expected for genuine conical 
waves 13]. Here we concentrate our attention in the 
beam-periphery structure in order to demonstrate that 
filaments are indeed conical and, more precisely, NL- 
UBB waves. To this end we modified the diagnostic 
by adopting professional digital photo-camera (Canon- 
EOS D30), which has the unique advantage of permitting 
strong local saturation (by the central spike) without any 
blooming effect. The experiment was done by launching 
a spatially-filtered, collimated, 200 fs, 0.1 mm, 1.5 /xJ, 
527 nm Gaussian wave packet into a 31 mm water-filled 
cuvette. By using an imaging spectrograph we verified 
that, at the specified pump energy, no relevant spectral 
broadening occurs, the self-phase modulation occurring 
in the sole spatial domain. Figure [3 shows the measured 
fluence profile at the output facet of the non-linear sam- 
ple, and a fitted UBB profile. BB and Airy profiles (only 
one half, for clarity) are also shown for comparison. One 
can appreciate: (i) The Bessel-like decay (~ 1/r) of the 
measured profile in a fairly vast region of the beam, which 
distinguishes it sharply from any (Airy-type) aperture- 



diffraction pattern (with faster decay ~ 1/f 2 )- (h) The 
accurate fitting of the UBB radial modulations to those 
of the measured profile in the full recorded area, which 
even reproduces the increasing frequency of the modu- 
lations (respect to a BB) towards the beam center, at- 
tributable to Kerr self-focusing, (iii) The reduction in 
modulation contrast (compared to a BB), which is a sig- 
nature of unbalance between inward and outward conical 
power flows. 

In conclusion, we have reported on the existence, char- 
acteristics, stability and experimental relevance of non- 
linear unbalanced-Bessel-beams (NL-UBB), i.e., the sta- 
tionary and localized solutions of the 2D+1 NLSE in the 
presence of NLL. NL-UBB are asymptotically linear coni- 
cal waves, carrying a net inward power flux that compen- 
sate for the NLL. We have shown that the UBB asymp- 
totics is the sole compatible with NLL, no matter which 
are the specific non-linear phase terms in the NLSE. 
Therefore NL-UBB are possible solutions also of nonin- 
tegrable NLSEs as well as of the Gross-Pitaevskii equa- 
tion for Bose-Einstein condensates. Owing to the spa- 
tiotemporal analogy, the results directly apply to pulse 
propagation in planar wave guides with anomalous dis- 
persion. More generally, similar waves should exist also 
in the case of 3D+1 NLSE with NLL, describing full spa- 
tiotemporal localization in bulk media with NLL. The 
unique property of long-range stationarity in the pres- 
ence of energy transfer to matter (or to other waves) 
makes NL-UBB ideal for several applications including 
deep-field non-linear microscopy, laser micro machining, 
laser- writing of channel wave guides, charged-particle ac- 
celeration, creation of long and stable plasma channels in 
atmosphere and, of course, energy transfer between dif- 
ferent types of waves (e.g., optical — > X, or optical — > 
teraHerz). 
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